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Abstract. This proceeding briefly summarizes our recent VISHNU hybrid model investigations 
on the chemical and thermal freeze-out of various hadrons species in 2.76 A TeV Pb+Pb 
collisions. Detailed analysis on the evolution of particle yields and the last elastic collisions 
distributions during the hadronic evolution reveals that the two multi-strange hadrons, H and 
fl, experience early chemical and thermal freeze-out when compared with other hadron species. 


1. Introduction 

In relativistic heavy-ion collisions at top RHIC and the LHC energies, a hot and dense matter 
- the Quark-Gluon Plasma (QGP) has been created. During its expansion, the QGP fireball 
quickly cools down, which then undergos the phase transition and produces a large amount of 
hadrons. These hadrons subsequently experience frequent inelastic and elastic collisions during 
the hadronic evolution. With the termination of inelastic and elastic collisions, the evolving 
system reach the chemical and thermal freeze-out, respectively. 

Traditionally, the chemical freeze-out temperature Tch and the baryon chemical potential 
are extracted from the particle yields of various hadrons using the statistical model mmm, which 
gives a uniform Tch and ^b for all hadrons species. Similarly, pure hydrodynamic simulations 
implement a uniform thermal freeze-out temperature to construct the freeze-out hyper-surface, 
which decouples various hadrons from the bulk matter [3]. 

Microscopically, the related inelastic and elastic scattering channels are different for each 
hadron species. As a result, the chemical and thermal freeze-out procedures are very possibly 
hadron-species dependent. In this proceeding, we will briefly summarize our recent investigation 
on the chemical and thermal freeze-out of various hadrons in 2.76 A TeV Pb-|-Pb collisions [S], 
based on dynamical simulations from VISHNU hybrid model. 

2. Setup of the calculations 

Our calculations implement VISHNU hybrid model [6] that combines (2-|-l)-d relativistic viscous 
hydrodynamics (VISH2+1) [7] for the QGP fluid expansion with a microscopic hadronic transport 
model (UrQMD) [8] for the hadron resonance gas evolution. The hydrodynamic calculations input 
the equation of state (EoS) s95p-PGE |9] with the transition to UrQMD at a constant temperature 
of 165 MeV. Eollowing Ref. [101 lllj . we input smooth initial conditions generated from the MG- 
KLN model [121 fT3] , where the normalization factor of the initial entropy density profiles is tuned 


Figure 1. (Color online)Particle 
yields of tt, K, p, A, H, and 
n in the most central Pb+Pb 
collisions at y/sj^=2.76 TeV. 
The experimental data are from 
the ALICE Collaboration [161 
[m HE]. Theoretical curves 
are from VISHNU hybrid model 
calculations, using MC-KLN initial 
conditions, rj/s = 0.16 and = 
165 MeV. 
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to reproduce the multiplicity of all charged hadrons in the most central Pb+Pb collisions |14j . 
The hydrodynamic simulations start at tq = 0.9 fm/c with the specific shear viscosity {p/ s)qgp 
set to 0.16, and the bulk viscosity set to zero [iniiis]. 

3. Results 

The chemical and thermal freeze-out of the evolving system are respectively related to the 
terminations of inelastic and elastic collisions. After the chemical freeze-out, the hadron yields 
no longer change. After the thermal freeze-out, the momentum distributions of various hadrons 
no longer change. In this section, we will briefly summarize our recent investigations on the 
chemical and thermal freeze-out of various hadron species at the LHC, based on dynamical 
simulations from VISHNU hybrid model. 

Figure 1 presents the VISHNU calculations of the particle yields for vr, AT, p, A, H, and H in the 
most central Pb+Pb collisions, which nicely fit the corresponding experimental data measured by 
ALICE. It was found that the baryon and anti-baryon {B-B) annihilations during the hadronic 
evolution reduce the proton and anti-proton {p/p) yields by ~ 30%, largely improving the 
description of the p/p data [H]. Ref. m and [5] also showed that VISHNU hybrid model also 
nicely fit the multiplicities, pp spectra and differential elliptic flow of vr, K, p, A, S, and H 
at various centrality bins. With the good descriptions of these related soft hadron data, it 
is meaningful to further investigate the chemical and thermal freeze-out procedures of various 
hadrons at the LHC. 

Figure shows the time evolution for the particle yields of vr, K, p, A, S and H during the 
hadronic evolution, simulated from VISHNU hybrid model with and without B-B annihilations. 
For the simulations without B-B annihilations, the yields of these two multi-strange baryons 
H and H almost do not change. We also notice that their yields quickly reach saturation for 
the case with B-B annihilations. This reveals that S and H experience early chemical freeze- 
out after the hadronazition. In contrast, the yields of p slowly decreases during the hadronic 
evolution especially for the case with B-B annihilations, indicating later chemical freeze-out 
of protons. On the production of A and K, the B-B annihilations almost balance with other 
inelastic collision channels. As a result, the yields of A and K only slightly change during the 
hadronic evolution with B-B annihilations [5]. 

The right panels of Fig. 1 presents the changing rates of the particle yields of these identified 
hadrons. For the simulations without B-B annihilations, the changing rates for K, p, and A 
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Figure 2. (Color online) Left 
panels (a)-(f): time evolution 
of relative particle yield density 
f (i)/f (0) for rr, K, p, A, =, 
and ri during the UrQMD expan¬ 
sion of VISHNU. Here, ^(t) and 

^(0) denote the particle yield 
density at mid-rapidity at later 
evolution time and at the start¬ 
ing time, respectively. Right pan¬ 
els (g) and (h): time evolution of 
the changing rates for the corre¬ 
sponding particle yield densities. 
Solid/dashed lines denote the 
VISHNU simulations with/without 
B-B annihilations. 


show wide peaks along the time axis, illustrating that the associated inelastic collisions are still 
frequent after 10-20 fm/c. Due to small hadronic cross sections, multi-strange hadrons H and D 
rarely interact with other hadrons. While, the annihilations with their own anti-particles lead to 
wide peaks of the changing rates for H and D in panel (h), which delay their chemical freeze-out, 
when compared with the case without B-B annihilations. However, these multi-strange hadrons 
still experience earlier chemical freeze-out in the full VISHNU simulations compared with other 
identified hadrons [5]. 

Figure [^presents the thermal freeze-out time distributions of vr, iF, p, A, H and D in the most 
central Pb-|-Pb collisions. Here, the thermal freeze-out time distributions are extracted through 
analyzing the space-time information of the last collisions for various hadron species. For H and 
D, the peaks of their thermal freeze-out time distributions are located around 10 fm/c. While 
the peaks of p and A distributions are around 20-30 fm/c. Due to a large amount of hadronic 
scatterings during the late UrQMD evolution, the thermal freeze-out time distributions of vr and 
K widely spread along the time axis. Fig. 3 thus illustrates that thermal freeze-out procedures 
are hadron species dependent. Compared with other hadrons, the two multi-strange hadrons H 
and D experience early thermal freeze-out due to their much smaller hadronic cross sections [5]. 

4. Summary 

In this proceeding, we briefly review our recent investigations on the chemical and thermal 
freeze-out of various hadron species in 2.76 A TeV Pb-|-Pb collisions within the framework 
of VISHNU hybrid model. We found that the chemical and thermal freeze-out procedures are 
hadron-species dependent. Due to their much smaller hadronic cross sections, the two multi- 
strange hadrons H and D experience earlier chemical and thermal freeze-out than other hadrons. 
To extract the possible effective chemical freeze-out temperatures of various hadron species 
requires detailed analysis of the space-time distributions of the last inelastic collisions, which we 
would like to leave it to future study. 




























Figure 3. (Color online) 
Thermal freeze-out time dis¬ 
tributions for vr, K, p, A, H 
and O in the most central 
Pb-|-Pb collisions, calculated 
from VISHNU with B-B anni¬ 
hilations. 
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